Effect of phonon scattering by surface roughness on the universal thermal conductance 
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The effect of phonon scattering by surface roughness on the thermal conductance in mesoscopic 
systems at low temperatures is calculated using full elasticity theory. The low frequency behavior 
' of the scattering shows novel power law dependences arising from the unusual properties of the 

elastic modes. This leads to new predictions for the low temperature depression of the thermal 
conductance below the ideal universal value. Comparison with the data of Schwab et al. [Nature 
404, 974 (2000)] suggests that surface roughness on a scale of the width of the thermal pathway is 
important in the experiment. 
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Thermal transport in mesoscopic systems at low temperatures shows universal properties analogous to the quantized 
electrical conductance jlj . Theoretical analyses of phonon transport when the thermal wavelength becomes comparable 
[ to the dimensions of the thermal pathway predict a thermal conductance K that is independent of many of the details 
jjy • of the geometry and material properties ||, and a universal value for K/T = NoTT 2 k^/3h at low enough temperatures 
H Q. (Here Nq is the number of modes with zero frequency at long wavelengths, equal to 4 for a single, free-standing 
elastic beam.) These predictions have since been connected to more general results on bounds on entropy transport at 
| low temperatures JsJ , and to thermal transport by particles of arbitrary statistics || . The recent confirmation of the 
3 ■ universal thermal conductance in tiny silicon nitride devices |?| is an experimental tour de force. Although verifying 
the predictions of a universal value of K/T at low enough temperatures, the experiments showed values of K/T that 
decrease as temperature increases in the range of 0.084K to 0.2K, before beginning to rise at higher temperatures as 
more vibrational modes that can carry the heat are excited. 

We H and others |9], [l0| have previously studied a simplified treatment of this problem using a scalar model for 
the elastic waves. However elastic waves in confined geometries have many unusual features, such as modes with 
a quadratic dispersion relation u> oc k 2 at small wave numbers k and regions of anomalous dispersion duj/dk < 0, 
that are not captured in this simple model and might be expected to have a strong influence on the low temperature 
J> ■ transport. 

In this letter we study the general problem of the scattering of phonons by rough surfaces for a rectangular 
• cross-section elastic beam using three dimensional elasticity theory. The scattering of the low frequency modes that 
ly— j \ contribute to the conductance at low temperatures depends on the detailed properties of the elastic modes, and we 
find novel power law dependences for the frequency dependence of the scattering off unstructured roughness that 
t— I | are not those anticipated by simple analogy with Rayleigh scattering. In turn this low frequency behavior yields a 
depression of K/T at low temperatures with unexpected power laws, with a faster decrease of K/T as the temperature 
is raised than anticipated from the scalar model. Comparing these predictions with the experiments of Schwab et al. 
suggests that unstructured roughness is inadequate to explain the data: instead we find a roughness length scale 
G \ comparable to the width yields a fit consistent with the experimental trends. 
*^ . In the ballistic transport regime, the thermal conductance at temperature T takes the form j|, H @ 

o 

> 
X 



where the sum is over the modes m propagating in the structure and oj m is the cutoff frequency of the rath mode, 
and (3 — 1/(&bT). The parameter T m (u)) is the transmission coefficient: T m (u) ^ 1 corresponds to a reduction in the 
transport due to some scattering process. At low temperatures, only the modes with uj m = contribute, and in the 
absence of scattering the conductance takes on its universal value K u . 

To calculate the transmission coefficient T m {ui) in Eq. ([!]) we use a Green function method similar to the one 
previously reported in the scalar- wave calculation || . The thermal pathway is modelled as a rectangular elastic beam 
or waveguide of width W, depth d, and length L in the x-direction. We assume smooth top and bottom surfaces 
corresponding to the epitaxial growth planes, and rough surfaces at y = ±W/2 + f±(x) deriving from the lithographic 
processing. The functions f±(x) define the roughness, which we assume to be independent of the z coordinate and 
small f± < W. 
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Assuming isotropic elasticity theory, the displacement field u satisfies the wave equation: 

pd 2 u t /dt 2 - dTij/dXj = 0, 



(2) 



where p is the mass density, and is the stress tensor, = CijkidkUi with Cijki the elastic modulus tensor. 
Stress free conditions apply at the surfaces T^n^s — with n the surface normal. The Green function Gij{x 1 x' ;t) 
is defined to satisfy the corresponding equation with a source term £jj(5(x — x')£(t — t') on the right hand side, and 
stress free boundary conditions at the smoothed surfaces y = ±W/2. Using the completeness of the orthonormal, 
smooth-surface, elastic modes ul™' (x, u>) at frequency u> (with k the wave number along the wave guide and n labelling 
the transverse mode structure), the Fourier transformed Green function can be written in the form 



(n)* I _./ \ (ri)/-> \ 
U] [X ,0J)U) {X,0J) 



o (™) 

ZpUJVg 



(3) 



The sum is over the elastic modes n that propagate at frequency uj: the condition u){k) = lj with to(k) the dispersion 
relation defines a discrete set of wave numbers k n and corresponding mode indices n. For monotonic uj(k) there is 
a single solution for k on each branch (defined by a continuous uj(k)), but with anomalous dispersion there may be 

(n) 

multiple solutions. v g is the group velocity of mode n at frequency ui. 

The stress free boundary conditions for the propagating waves apply at the actual (rough) surfaces; expanding in 
small f± lead to effective stress boundary conditions at the smooth surfaces y — ±W/2, which can then be included 
using Green's theorem. These manipulations lead to the expression for the scattered field u' sc ' to lowest order in 
the small roughness in terms of the stress field of the incident wave at y = ±W/2. The back-scattered field from a 
unit-amplitude incident wave in mode m is 
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with 
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(repeated indices i and ± are to be summed over) where E is Young's modulus, p — E/2 (1 + a) the shear modulus, 
a the Poisson ratio, and is the stress field corresponding to the displacement field u^ n -*. The mode sum in Eq. (jij) 

is confined to modes propagating to negative x, i.e. Vg (k n ) < 0: usually these will be modes with negative wave 
numbers k n , but there are also regions of anomalous dispersion which have a negative group velocity for positive k n - 
Notice that Eq. (Q) is comprised of separate kinetic, compression and shear terms. 

Although the expression Eq. (|^) applies quite generally to a three dimensional waveguide with a rectangular cross- 
section [p~5[| , there are no closed- form expressions for the modes in general, and so a direct evaluation of the 
scattering would have to be done with numerical solutions for the modes. Here, we instead use the thin plate 
approximation d <C W [p"2[ |l3f| , which yields closed form expressions for the displacement fields of the modes (in terms 
of a dispersion curve u>(k) given by numerical solution of a simple transcendental equation). This approach correctly 
accounts for the important properties of the elastic modes, for example exactly reproducing the small k dispersion 
relations including the quadratic dispersion of the bend modes, and showing regions of anomalous dispersion, whilst at 
the same time allowing analytic progress . Thin plate theory should be quantitatively accurate for many mesoscopic 
experiments where the thickness of the sample (formed by epitaxial growth) is often much less than the width (given 
by lithography). The dispersion curves in the thin plate approximation are shown in Fig. (Q). Note that the modes 
polarized in the plane and perpendicular to the plane (flexural modes) do not couple. 

An important simplification in thin plate theory is that all z-components of the stress tensor Tj Z , T Z i in Eq. (^) may 
be put to zero 12 . The expression for the flux scattering rate (per unit length of rough waveguide) from mode m to 

mode n at frequency u>, defined as v g n ^ /v g m ^ times the ratio of the mode intensities, simplifies to 
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FIG. 1: Dispersion relation for in-plane modes (solid) and flexural modes (dashed) for a geometry ratio d/W = 0.375 and 
Poisson ration 0.24. The wave numbers are scaled with the width W, and the frequencies by W/ct with ct = y /J./ p. 
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TABLE I: Scattering coefficients for the zero onset frequency modes at low frequencies: c denotes compression, b denotes bend, 
t denotes torsion, bb denotes bend to bend scattering etc. Values are quoted for ■y m W 4 /g(0) as a function of scaled frequency 
uj = loce/W . For the flexural bend to bend scattering (bb) the terms in the braces in Eq. (H) cancel to leading order resulting 
in very small 0(lj 3 ) scattering. There is no scattering between in-plane and flexural modes for the z-independent roughness 
assumed. 

where 4> gives the transverse dependence of the modes it^™' = e lknX 4>[ n \y, z), and again back-scattering corresponds 
to modes n with negative group velocities which except for the regions of anomalous dispersion correspond to negative 
k n . To arrive at this expression we have performed an ensemble average over the surface roughness function f±(x), 
and g(k) is the spectral density, i.e. the Fourier transform of the roughness correlation function g{x) defined by 

(f(x)f(x'))=g(x-x'), (7) 

where / is /+ or /_, and we assume independent scattering off the two rough surfaces introducing a factor of 2 
into Eq. (j^). For weak scattering the transmission coefficient appearing in the equation for the thermal conductance 
Eq. (Q) can be obtained as T m = e~ lmL with j m = tI (») <0 7m,nj where we sum only over the back scattering, since 
scattering into forward propagating modes does change the heat transport. 

We first study the scattering behavior for unstructured "white noise" roughness, g(k) = g(0) independent of k. At 
low frequencies only the four lowest modes with zero onset frequency survive (the compression mode and a bending 
mode polarized in the plane, and the torsion mode and a second bending mode polarized perpendicular to the plane). 
Explicit expressions for the scattering coefficients from these modes 7 rn ,(w)W /4 /<7(0) are shown in Table [j| as a function 
of the scaled frequency u> = loce/W with ce — \J E/p the propagation speed of extension waves in the beam. Notice 
the back-scattering from the compression mode to the time reversed mode, and from the torsion mode to its time 
reversed mode, have the uj 2 dependence anticipated in analogy with Id Rayleigh scattering (as indeed was found in 
the scalar- wave calculation). On the other hand scattering to or from the bending modes introduces quite different 
power law dependences, that will dominate at low frequencies. Most of this novel behavior can be understood simply 
from the prefactor in Eq. (0), 7 oc u> 2 / Vg™ 1 ^ Vg n \ since v g goes to a constant at small frequencies for the compression 
and torsion modes, and is proportional to lu 1 ^ 2 for the bending modes. The remaining part of Eq. (^), the integral 
over z, is 0(1) as uj — > 0, except for the flexural bending mode back-scattering where the two contributions cancel at 
leading order. Note also that the expressions for the flexural modes involve additional factors of the width to depth 
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ratio W/d, so that in the thin plate limit these modes will be scattered more strongly at a given ui. 

The low temperature thermal conductance can be derived directly from the low frequency scattering expressions: 
if we write these as 7L = AQ P ', then the corresponding contribution of the suppression of the thermal conductance is 

5K/K U = AI p (T/T e )p (8) 

with Te = hcE/ksW and I p the integral 

3 r°° y p+* e y 

giving the same power law for the temperature dependence. 




scaled frequency, roW/ c t 

FIG. 2: Scattering coefficient 7 m W /4 /p(0) of the lowest (zero onset frequency) modes to all other modes for unstructured 
roughness: dashed-dotted - compression; solid - in-plane bend; dashed - flexural bend; dotted - torsion. 



Figure (|2|) shows the scattering coefficients over an extended frequency range, again with g(k) = const. The phonons 
are strongly scattered at each mode onset frequency, as a result of the zero group velocity, and the correspondingly 
large density of states to scatter into. We find additional strong scattering in the regions of anomalous dispersion 
(e.g. ujW/ct — 5), since again these portions of the dispersion curve are quite flat. Also, since the group velocity 
of the in-plane bending mode approaches zero, the scattering of this mode increases relative to others as u> — > 0. 
This agrees with the low frequency analysis, 7 oc u>. (The back-scattering of the flexural bending mode is reduced 
by the cancellation between the kinetic and stress scattering.) For intermediate frequencies 0.08 < uiW/ct < 1-8 the 
scattering of the torsion mode is larger than for the other modes. 

With this detailed understanding, we now attempt to fit the experiments of Schwab et al. We focus on the data 
below 0.4K where thin plate theory captures the relevant modes ]Tq| . To compare with experiment we first subtract 
the measured conductance from the ideal conductance predicted using the cutoff frequencies calculated numerically 
using the xyz approach (ll[] using d = 60nm and W = 1 60nm: this gives us an estimate of the effect of scattering 
which we can compare with our calculated scattering. The data processed in this way, Fig. (^|), shows an abrupt 
onset of scattering at some non-zero temperature (T ~ 0.08K): this does not appear to be consistent with the power 
law behavior predicted using unstructured roughness, Table | and Eq. (||). Notice that the discrepancy is worse than 
anticipated in the scalar-wave calculation Q, which gave larger exponents for the power law dependence. The delay 
in the onset of scattering suggests that the phonons with long wavelength are not much affected by the rough surfaces. 
To accommodate this fact, we have investigated the scattering by roughness centered around a finite length scale k^ 1 , 
using the parameterization 

g(k) = V^5 2 ae- a2 ( fc - fc °) 2 (10) 

For koa > I this results in strongly reduced scattering at long wavelengths. There are three roughness parameters 
that now characterize the surfaces: the roughness amplitude S, the correlation length a, and shift Ajq. We obtain a 
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FIG. 3: Thermal conductance change relative to the universal value K u as a function of temperature: circles - data of Schwab et 
al. subtracted from the predicted ideal conductance; solid line - predictions from the elasticity calculation. The insert shows the 
thermal conductance relative to the universal value K u as a function of temperature: circles - data of Schwab et al.; dotted line 
- ideal conductance calculated using the xyz method; and solid line - fit from elasticity theory with scattering. The roughness 
parameters used were 5/W = 0.11 and koW = 4.7. 



reasonable fit to the low temperature onset of the scattering using parameters a/W = 3, 5/W — 0.11, and k$W = 4.7. 
These parameters correspond to significant roughness at length scales comparable to the width of the device, which 
appears consistent with the electron- micrograph of the structure used [Q . 

In summary, we have examined the scattering of phonons by surface roughness and the effect on the universal 
thermal conductance. At low temperature, the scattering shows strong dependence on mode structure. In this 
temperature range, the elasticity model provides a better understanding of the scattering behavior and a reasonable 
fit to the experimental data if we assume that the roughness is concentrated at a length scale comparable to the 
lithographic scale. Further experiments at lower temperatures would provide a useful test of the predictions of the 
theory. 
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